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Abstract 

We present a theoretical study of transport properties of a liquid comprised of particles inter- 
acting via Gaussian Core pair potential. Shear viscosity and self-diffusion coefficient are computed 
on the basis of the mode-coupling theory, with required structural input obtained from integral 
equation theory. Both self-diffusion coefficient and viscosity display anomalous density dependence, 
with diffusivity increasing and viscosity decreasing with density within a particular density range 
along several isotherms below a certain temperature. Our theoretical results for both transport 
coefficients are in good agreement with the simulation data. 
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I. INTRODUCTION 



Liquids comprised of particles interacting via a Gaussian Core (GC) model pair po- 
tential have been receiving a lot of attention recently, [l- 16 1 Strong theoretical interest in 
GC model stems from the fact that this bounded potential is useful for describing interac- 



tions between inherently penetrable entities, such as polymer coils. |l3l . As such, both 



t hermo dynamic 



jsl, [igI and dynamic js-llll properties of GC fluids and binary mixtures have 



been extensively studied. It has been found that their transport coefficients exhibit anoma- 



lous behavior strongly reminiscent of waterlike model systems |17H2q. with diffusivity in- 
creasing and viscosity decreasing with density over a certain range of thermodynamic con- 
ditions. Furthermore, a strong correlation between this behavior and structural anomalies 
quantified via excess entropy has been demonstrated, jo- 11 1 Both structural and dynamical 
anomalous behavior has been rationalized in terms of bounded nature of GC potential which 
results in increasing amount of interparticle overlap with increasing density. In view of the 
above, it would be of interest to provide a firm theoretical link between structural and trans- 
port properties of GC fiuid on the basis of microscopic statistical mechanical theory. Such a 
connection has been recently established for waterlike model systems j^] by combining inte- 
gral equation theory of structure with mode-coupling theory (MCT) treatment of dynamics. 
The goal of the present work is to develop a similar treatment for the GC model. We show 
that integral equation theory/MCT combination is indeed capable of capturing anomalous 
behavior of transport coefficients of both neat GC fiuids and binary mixtures observed in 



molecular dynamics (MD) simulations, fl-lnl Our MCT-based microscopic analysis helps to 



shed further light on the origin and nature of transport anomalies. In particular, we are 
able to rationalize why viscosity anomaly persists over much more narrow range of densities 
and temperatures compared to the diffusion anomaly. 

The remainder of the paper is organized as follows. In Section |IT1 we describe the micro- 
scopic model and review the MCT approach employed to calculate the transport coefficients. 
In Section Hm we compare theoretical results with the simulation data for the shear viscosity 
and the self-diffusion coefficient. In Section [IV] we conclude. 
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II. MICROSCOPIC MODEL AND THEORY 



We consider a system comprised of spherical particles interacting via isotropic GC pair 
potential 0(r): 

0(r) = eexp [-(r/fj)^] , (1) 

where the two parameters e and a characterize the height and the width of the interaction 
profile, respectively. In the previous studies of the GC model, integral equation theory 
with the hypernetted chain (HNC) closure was shown to give reliable structural results for 
a wide range of densities studied, Hence, we employ the HNC closure throughout this 
study to compute the radial distribution function g{r) for a liquid described by the GC 
pair potential. The validity of this approach will be further confirmed by comparing various 
structural quantities with the corresponding simulation data. 

The main focus of the present work is the calculation of the transport properties of GC 
fluid and analysis of their anomalous behavior. We first describe our treatment of the self- 
diffusion coefficient. The latter is obtained from the total time integral of the time- dependent 
friction C(^)- 



27 



= ^. (2) 
with 

Co= / dtat). (3) 
Jo 

The MCT result for the time-dependent friction reads: 

C W = ^ r dkk'c{kfF{k, t)Fs{k, t) (4) 

where m is the mass of the fluid particle, T is the temperature, p is the number density, and 
c{k) is the fluid direct correlation function, which we obtain from the HNC closure. In the 
above, F{k,t) is the fluid dynamic structure factor, which we compute from the continued 



fraction representation of its Laplace transform truncated at the second order: |27l. l29| 



z+T-\ky 

where S{k) is the fluid static structure factor, and 5i{k) is the initial time value of the 
order memory function (MF) of F{k, t). For the parameter T~^{k) we use the expression due 



to Lovesey: [sol r ^(fc) = 2-^52{k) /t: . The quantities 6i{k) and 52{k) can be easily calculated 
from the first three short-time expansion coefficients of F{k,t); the microscopic expressions 
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3l| Finally, Fs{k,t) 



for the latter are well-known and will not be reproduced here, 
is the fluid self-dynamic structure factor, for which we have adopted a simple Gaussian 
modeh[27|,|28| 

-Ko _ V 



Fs{k,t) = exp 



-ksTk^ ( 1 



mCo 



Co 



(6) 



Given that the self-dynamic structure factor is a function of Co? which, in turn, depends 
on Fs{k,t) via Eq. (jl]), Eqs. (I4])-(I6]) need to be solved iteratively and self-consistently. One 
could use a more accurate model for Fs{k, t) in terms of the velocity time correlation function 
(TCF) of a tagged fluid particle, 321] but our numerical calculations have shown that this 
does not change the results for D in a noticeable way. 

We note that the MOT result for the time-dependent friction given by Eq. (jlj) arises from 
coupling to collective density modes, and, as such, is only expected to be valid at intermediate 
to long times, i.e. it describes the slowly varying tail of C{t). At short times, one typically 
adopts an empirical approach by modeling the initial decay of the time- dependent friction 
Via so,„e ..ap,d.y decaying analyfca. f„,>ct,o,> (e.g. a Ga„ss,apl which ,s constructed in sud, 



a way as to preserve the exact short-time 
result for the initial time value C(0) reads: 

Airp 



Dehavior of C(^)- [27|, |29|] In particular, the exact 



27, 



29] 



C(o) 



3m 



drr'^g{r)'V'^(f){r). 



(7) 



At the same time, the zero-time value of the MOT expression for the time- dependent friction 
is given by: 



C(o)^ 



ksTp 



dkk^c{kfS{k) 



(8) 



By comparing the two expressions given by Eqs. (JTj) and ([8]) and by applying Parseval's 
theorem, one can see that the two results are equivalent under the random phase approx- 
imation (RPA), c(r) = —(j){r)/kBT. It has been shown in the previous studies of the GC 
liquid[l3| that RPA is quite accurate in describing its structural properties. As such, MCT 
method provides essentially correct zero-time value of the time- dependent friction. Further- 
more, numerical results presented in the next Section demonstrate that the MCT approach 
also describes the initial decay of C,{t) quite accurately. Hence, in the present work we do 
not decompose time-dependent friction into binary and collective terms, but rather employ 
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Eq. dl]) at all times. As will be seen below, a similar situation occurs in the MCT treatment 
of the potential part of the shear stress autocorrelation function (SACF), which enters into 
the calculation of the shear viscosity. 

The microscopic expression for the shear viscosity is given by the Green-Kubo formula in 
terms of the total time integral of the SACF (i.e. time- dependent shear viscosity r7(t)):[33l] 

V = TTT^ / dt{J^y{t)J^y{0)) = / dtr]{t). (9) 



VkeT 

with the off-diagonal components of the stress tensor J^.^ given by: 

N N N ^ . . 

J.y = E mv^vf -J2Y: -l-^ = J'v + J'^y (10) 

i=l i=l j>i ^1 

In the above, V is the total volume, N is the number of particles, lower indices i and j 
label the particles, while upper indices x and y denote the vector components of the particle 
velocities and displacement vector r^j connecting the particles i and j. 

According to Eq. ( ITOl) . the off-diagonal elements of the stress tensor are comprised of 
kinetic and potential terms, J^y and J^y, respectively. As a result, the SACF splits into 
three individual parts: potential-potential, kinetic-kinetic, and mixed kinetic-potential con- 
tributions: 

ri{t) = Vppit)+Vkk{t)+2vkp{t) = ^^{{Jly{t)Jl^m + (4(t)J4(0)) + 2(4(t)J4(0))} , 

(11) 

with the shear viscosity coefficient given by the sum of the total time integrals of the three 
terms above: 

Tj = rjp + rjk + 2r]kp, where rjp = dtr]pp{t) etc. We now discuss the calculation of each of 
these terms in turn. 

As in the case of the time-dependent friction, the MCT treatment of the potential term 
rjfmit), involves coupling to collective density modes, and the corresponding result reads: 



29| 



Sik) 



2 



(12) 



where prime denotes differentiation with respect to the argument. In general, the above 
expression is used at intermediate and long times, while short-time behavior of rjppit) is 
modeled phenomenologically, based on its short-time expansion coefficients. The exact result 
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for the initial time value of ripp{t) is given by:|27l. |29[ | 

^..(0) = -p drgir)- 



dr 



(13) 



while the corresponding MCT result can be obtained from Eq. ( !T2|) : 



dkk^ 



S'{k) 



[Sik) 



(14) 



GOtt^ jo 

By comparing the two expressions given by Eqs. ( IT3l) and (fT4l) and by applying Parseval's 
theorem, one can see that the two results agree when RPA holds. As will be seen from our 
numerical results presented in the next Section, not only zero-time value of rjppit), but also 
its initial rapid decay is well reproduced by the MCT. Hence, we use Eq. (fT2l) to model the 
potential part of the time-dependent shear viscosity at all times. 



Next, the MCT expression for the mixed kinetic-potential term is given by: 34l | 



Vkp{ty 



m 



dF{k,t) 
di 



(15) 



IStt^ Jo ' S{ky 

In this case, the initial-time value %p(0) = 0, while the first (linear) term in the short- 
time expansion of i]kp{t)^^'^ agrees (under the RPA) with the corresponding exact expansion 
coefficient. 35|] Once again, we employ the MCT expression to describe fjkpit) at all times. 



Final 
reads: [34 1 



y, the MCT result for the kinetic contribution to the time-dependent friction 



m 



dkk^ 



-Cl{k,t) + \cl{k,t) + Cu{k,t)Cn{k,t) 

6 



(16) 



The upper cutoff on the wavevector integral, k^ax = (Gvr^p)^/'^, is determined by the 
requirement [34 ] that the initial time value rjkkiO)^'^^ coincides with the exact value rjkkiO) = 
pksT. In the above, Cii{k,t) = —F{k,t)/k'^ is the longitudinal current TCP, and Ctt{k,t) is 
the transverse current TCP. The Laplace transform of the latter (truncated at the second 
order) is given by: [271, l29| 

Ctt{k,Z) = ^— ryr (17) 



z + 



z+ 



z+T,-\ky 

where 5it{k) is the initial time value of the order MP of Ctt{k,t). The quantities Sit{k) 
and S2t{k) can be obtained from the first three short-time expansion coefficients of Ctt{k,t), 



which are well-known. 



In analogy to the continued fraction representation of F{k, t), parameter ^{k) is taken 



to be proportional to a/ S2t{k): T^ ^{k) = ^^/S2t{k), with the proportionality factor £ deter- 
mined by requiring that Ctt{k,t) behaves correctly in the hydrodynamic limit: 271. |29 



[Ctik.t)] 



knT 



hyd 



exp 



k'^r]t 
pm 



(18) 



where rj is the shear viscosity coefficient. The above requirement yields: 



(19) 



pm fe-)-o 6it{k) 

Thus, the calculation of the transverse current TCF requires the knowledge of the shear 
viscosity coefficient (which is needed to compute parameter ^ above), while t] itself depends 
on Ctt{k, t) via the kinetic part of the time- dependent shear viscosity. Hence, the calculation 
of the transverse current TCF and the viscosity coefficient needs to be performed iteratively, 
in analogy to the calculation of the self- diffusion coefficient and time-dependent friction 
described earlier. 

In order to assess the accuracy of the MCT approach described above, in the next section 
we compare our theoretical results with the MD datajs-?, 9-ll| for self-diffusion and shear 
viscosity coefficients of GC fluid. 



III. RESULTS 



Our results will be presented in terms of dimensionless density and temperature deflned 
by: p* = pcr^ and T* = ksT/e. We also deflne dimensionless variables for time, self-diffusion 
coefficient, and shear viscosity coefficient as follows: t* = t{e/ma'^Y^'^, D* = D(m/ea'^y^'^, 
and rj* = rja^/ ^me. 

Before presenting and discussing our results for transport coefficients, we ascertain the 
accuracy of the HNC closure in calculating structural properties of GC fluid. This is of par- 
ticular importance because earlier studies have revealed a deepconnection between transport 
coefficients and structural properties, such as excess entropy. [9|4ll| To this end, we compute 
from integral equation theory with HNC closure gir) for the GC fluid, from which one can 
readily obtain the two-body contribution to the excess entropy deflned as follows: 



§2 = — 27rp / drr^ [5'('") In (?(?") — (fi'('") — 1)] 
'o 



(20) 
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MD dataj^ and our HNC results for —S2 are presented in the upper panel of Fig. [T] as a 
function of fluid density along two isotherms: T* = 0.08 and T* = 0.2. In general, theory is 
in good agreement with simulations except in the high density range, where it overestimates 
—S2 somewhat. Both MD and HNC results show that the negative two-body excess entropy 
initially increases with density (up to p* ~ 0.25), reflecting increasing structural order of 
the fluid upon its compression. Concomitantly, radial distribution function (not shown) 
progressively becomes more structured, with the height of its first peak increasing and the 
second solvation shell gradually developing. 5, ol In simple fluids, whose interaction potential 
contains a steeply repulsive core, this familiar type of behavior typically persists throughout 
the entire liquid density range. By contrast, GC fluid displays qualitatively different behavior 
at higher densities, whereby its structural order (as measured by — S2) decreases upon further 
compression, which manifests itself in flattening of solvation shells in 5'(r).j^, 0] This is a 
consequence of the bounded nature of the GC potential, which allows interparticle overlap. 
The latter becomes more prominent at higher densities and/or temperatures, ultimately 
resulting in a high density ideal gas-like structure, [l, 3| 

In the high-density regime, GC system behaves as a weakly correlated "mean field 



fluid" 



13| whose direct correlation function is adequately described by RPA approach. This 



fact can be illustrated by comparing the results for the initial time value of the time- 
dependent friction given by Eqs. ([7]) (exact) and ([8]) (MCT), which would yield identical 
results if RPA were exact. The results from the two methods for the dimensionless zero- 
time value C*(0) = C(0)mcrVe are shown in the middle panel of Fig. [1] as a function of fluid 
density along the isotherm T* = 0.08. Indeed, the two approaches give quite similar values, 
and the difference between them diminishes with increasing density. In addition, we have 
applied the same two methods to compute the short-time expansion coefficient of the 
(normalized) time dependent friction, C2(0) = —^-^^\t=o/({^)- The corresponding dimen- 
sionless result, Qi^) = C2(0)mcr^/e, is presented in the lower panel of Fig. [H One sees that 
the two values are remarkably close numerically (except in the low density regime), even 
though the two expressions (exact and MCT) do not become equivalent under the RPA in 
this case. 

By comparing the density behavior of ({0) and C2(0) shown in Fig. [H one observes that 
the latter increases nearly linearly with density throughout the entire range, while for the 
former fast increase at low densities is followed by slower growth at intermediate p and 
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Figure 1 



1 1 ' ' 

• T =0.08 (MD) 

- T*=0.08 (HNC) 

^ T*=0.2(MD) 




FIG. 1: Upper panel: Simulation and theoretical results for the negative two-body contribution to 
the excess entropy as a function of fluid density along two isotherms. Middle panel: Initial value of 
the time-dependent friction as given by exact expression (Eq. d?])) and the MCT approach (Eq. ([8])). 
Lower panel: Short-time expansion coefficient of order t"^ of the normalized time-dependent friction 
from exact and MCT approaches. 

nearly flat behavior in the high-density regime. These results signify that the initial decay 
rate of the time- dependent friction grows continuously with density, while the growth of its 
zero-time value with p gradually slows down and nearly saturates. 

The above observations can be used to rationalize the density behavior of C„(t), the 
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Figure 2 



T =0.08 




p =0.1 (MD) 
p*=0.1 (MCT) 
p*=0.3 (MD) 
p*=0.3 (MCT) 
p =0.3 (Ansatz) 
p*=2 (MD) 
p*=2 (MCT) 



FIG. 2: Simulation and theoretical results for the normalized velocity TCF of a tagged GC particle 
at T* = 0.08 at three different densities. 



normalized velocity TCF of a tagged GC particle defined by: 



m 



(21) 



where Vq is the x-component of the tagged particle velocity, and (■ ■ • ) denotes the Boltzmann 
equilibrium average. The Laplace transform of this TCF, Cy{z), can be related to the Laplace 



transform of the time dependent friction as follows: 

1 



27 



29| 



CJz) 



(22) 



Simulation [5|] and MCT results for Cy{t) at three different densities along the isotherm 
T* = 0.08 are shown in Fig. [2l The most notable feature is biphasic decay of the TCF: 
except at the lowest density shown, rapid initial decay of Cy{t) is followed by a "bump" at 
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intermediate times and a slowly decaying long-time tail, whose amplitude increases markedly 
with density. The initial decay of the correlation function is governed by its t"^ short-time 
expansion coefficient, or, equivalently, by the initial time value of its MF, i.e. time dependent 
friction. As can be seen from the middle panel of Fig. [H this value initially increases with 
density, and then gradually saturates. Accordingly, in going from p*=0.1 to p*=0.3 the 
short-time decay of Cy{t) becomes considerably faster, while for densities beyond p*=0.3 it 
stays essentially unchanged. At the same time, as follows from the lower panel of Fig. [T|, the 
short-time decay of ({t) itself becomes faster with p throughout the entire density range. 
Faster decaying MF corresponds to slower decaying TCF, which indeed manifests itself in 
the amplitude of the slowly decaying tail of Cv{t) growing with density. 

In comparing theoretical results with the simulation data presented in Fig. one sees 
that the MCT approach reproduces biphasic decay of the velocity TCF quite well, with the 
only significant discrepancy observed at intermediate times for p*=0.3. In addition to the 



MCT method, we have also explored an alternative approach to constructing C„(t).|35l. l37| 



This semi-phenomenological approach is based on the short-time expansion of the TCF and 



35, 3 



is frequently employed in computing transport coefficients via Green-Kubo relation. [35|, 
In particular, one empirical form that is frequently used to model the normalized TCF is 
written as follows: Cy{t) = cos(6t)/ cosh(at), where parameters a and b are chosen in order 



to reproduce the exact short-time behavior of Cy{t) up to the term of order t^. 35|, |37|] This 
ansatz has been successfully applied to calculate diffusion coefficient and other transport 
properties of simple atomic liquids, whose interaction potential contains a steeply repulsive 



term at short separations. 35|, 1371] However, inspection of Fig. [2] reveals that this particular 
form would not be appropriate in the present it would not be able to capture biphasic 

decay of the TCF. Instead, we have attempted to model the time-dependent friction via the 
above phenomenological form: 

, , cos(6t) 

m = C — ^ 23 
cosh(arj 

Such approach requires the knowledge of the short-time expansion coefficients of ({t) up to 
the order of t^, or, equivalently, the expansion coefficients of Cy{t) up to the order of t^. The 
latter are well-known and will not be reproduced here. 



For the density value of p*=0.3, we display in Fig. [21 alongside with the MCT result, 
Cv{t) obtained from the time- dependent friction constructed via Eq. (!23|) . One sees that this 
phenomenological approach agrees with the MD data even better than the MCT method, 
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both at short (by construction) and at intermediate times. The performance of ansatz at 
other densities is equally good and is not shown to avoid overcrowding the graph. 

As one last remark concerning velocity TCF, we note that its pronounced biphasic be- 
havior would seemingly sug gest applying MCT approach either to C^(t) itself, 39|] or to 



its second-order MF. 



401 . |4l| We have carried out such calculations, and found, somewhat 



surprisingly, that the results were less accurate compared to the approach based on the 
first-order MF outlined above. 

We next compute the self- diffusion coefficient of GC fluid via Eq. ([2]), with time- 
dependent friction obtained both from MCT (via Eq. (jl])) and from ansatz (via Eq. (1251) ). 
Our results for D* as a function of GC fluid density, together with MD data,j3, ^ are shown 
in Fig. ([3]). MCT is generally in good agreement with the simulation, except that it overesti- 
mates the value of D at intermediate densities at T* = 0.08, as one could already expect from 
comparing MCT and MD results for the velocity TCF presented in Fig. ([2]). By contrast, 
ansatz systematically under-predicts D at low and intermediate temperatures. Overall, both 
methods correctly reproduce the trends in the density dependence of the self-diffusion co- 
efficient observed in the simulations. Specifically, at low and intermediate temperatures, D 
initially decreases with p, passes through a minimum (at about the same density, for a given 
temperature, where —S2 passes through a maximum), and then keeps increasing throughout 
the remaining density range. For the highest isotherm shown (T*=1.5), the minimum in D 
is barely perceptible, and at still higher temperatures it disappears altogether (likewise for 
the maximum in — S2). As has been already pointed out, [9] a strong correlation between 
anomalous density dependencies of —S2 and D points to the structural origin of the transport 
anomalies of GC fluid. This conclusion is further re-inforced by the observation that the 
density anomaly in D is successfully captured by the ansatz expression for time dependent 
friction, which is constructed exclusively from the short-time expansion coefficients of C('^)- 
The latter are simply equilibrium averages of certain functions of the interaction potential 
and its derivatives, i.e. purely structural properties. Some dynamical information does enter 
the MCT approach, e.g. via the dynamic structure factor, but, again, F{k,t) was obtained 
from its continued fraction representation (Eq. ([5]), which is completely determined by the 
short-time expansion coefficients of dynamic structure factor. 

Gradual disappearance of the anomalous density behavior of self diffusion coefficient at 
higher temperatures can also be understood in terms of density dependence of short-time 
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FIG. 3: Simulation and theoretical results for the self-diffusion coefficient of a tagged GC particle 
as a function of fluid density along three isotherms. 

expansion coefficients. In particular, for T* >1.5, the initial value of time dependent friction 
no longer saturates at high densities, but keeps increasing throughout the entire density 
range (not shown). As a result, the total time integral of ({t) keeps increasing, meaning 
that D decreases monotonically with p. 

Returning to low and intermediate temperatures, the anomalous increase of self-diffusion 
coefficient with density (and decrease of structural order) is by no means unique to GC 



fluid; this behavior has been observed for several waterlike model systems, 17H26| as well 
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Figure 4 




* 
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FIG. 4: Simulation and theoretical results for the shear viscosity of GC fluid as a function of fluid 
density along three isotherms. 



421 ] An important difference 



as for colloidal systems with short-range attractive interactions, 
however, is that for these systems the anomalous structural and dynamical behavior is 
observed for a limited density range only, after which the system reverts to normal behavior. 
This is due to the fact that all these model potentials contain a short-range steeply repulsive 
term, which manifests itself in both structure and dynamics at sufficiently high densities. 
By contrast, GC potential is bounded, and for this model anomalous behavior persists till 
the highest density studied. 
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We next turn to the discussion of shear viscosity coefficient. Simulation {g, 7| and theo- 
retical results for the shear viscosity of GC fluid as a function of fluid density along three 
isotherms are shown in Fig. HJ Also presented are the MCT results for the components of rj 
arising from the potential term, rjp, and from the sum of kinetic and mixed kinetic-potential 
terms, rj^ + 2rikp. At the lowest temperature studied, T* = 0.01, the shear viscosity coeffi- 
cient displays a pronounced anomalous density behavior. Initially, t] increases rapidly with 
p, passes through a maximum around p* ~0.25, exhibits anomalous decrease with density 
until about p* ~1, and then reverts to normal behavior, i.e. increases with density. This 
is different from the diffusion anomaly, where D, after passing through a minimum in the 
low-density region, keeps increasing with p throughout the entire density range studied. In 
addition, the temperature range where viscosity anomaly is observed is substantially more 
narrow compared to the diffusion anomaly. Thus, along the isotherm T* = 0.025, the mini- 
mum and subsequent maximum in t] are barely perceptible, both in MD and MCT results, 
while at still higher temperatures shear viscosity coefficient increases monotonically with 
density. 

If one considers the decomposition of rj into its components, one observes that the anoma- 
lous behavior at the lowest temperature is dominated by the potential contribution. Fur- 
thermore, after passing through a maximum around p* ~0.25, rip keeps decreasing with 
density for all remaining values of p. By contrast, rik + 2r]kp term is monotonically increasing 
with density. While its magnitude is negligible at low and intermediate densities, it becomes 
comparable to rip in the high density region, which produces a minimum in total rj. As 
the temperature is increased, the relative importance of rjk + 2rikp term also increases, and 
eventually its contribution to r] becomes dominant, which results in a monotonic increase of 
the shear viscosity coefficient with density. 

Some further light on the origins of the anomalous density behavior of viscosity at low 
temperatures can be shed by considering the time dependent viscosity given by Eq. ( fTTI) . 
The corresponding simulation and theoretical results are shown in Fig. [9] for four different 
densities along the isotherm T*=0.01. One notices that in going from p* = 0.1 to p* = 0.3, 
the initial time value //(O) increases dramatically, while the decay rate of r]{t) grows only 
slightly. As a result, rj, which is given by the total time integral of ri{t), grows signiflcantly 
with density in this range. The situation is reversed in the intermediate density range, from 
p*=0.3 to p*=l. Here the initial time value grows weakly, while the decay of ri{t) becomes 
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Figure 5 




FIG. 5: Simulation and theoretical results for the time-dependent shear viscosity of GC fluid at 
r*=0.08 for four diff'erent densities. 

much faster, thereby producing anomalous density behavior of viscosity in this range. Fi- 
nally, for densities higher than p*=l, the initial time value keeps increasing gradually, the 
initial decay rate does not change significantly, and the amplitude of the long-time tail of 
r]{t) due to mode-coupling effect grows steadily. Hence, the total time integral of ri{t) grows 
with p, and the shear viscosity coefficient reverts to normal density behavior. In terms 
of comparison between theory and simulation, we note that MCT underestimates the am- 
plitude of the long-time tail somewhat, which results in the under-prediction of t] seen in 
Fig.H 

In addition to the density behavior of shear viscosity at a given temperature, it is also of 
interest to analyze its temperature behavior at a given density. To this end, we present in 
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Figure 6 




FIG. 6: Simulation and theoretical results for the shear viscosity of GC fluid as a function of fluid 
temperature along three isochores. 

Fig. [6] simulation {g, 7| and theoretical results for the shear viscosity of GC fluid as a function 
of fluid temperature along three isochores. For the isochore p*=0.3 (where a maximum in rj 
as a function of p is observed), viscosity decreases with temperature at low and intermediate 
densities. The initial sharp decrease in this T range is associated with moving away from 
the liquid-solid phase boundary, js, 7| and the behavior of r] is dominated by the potential 
term. At higher temperatures, r]k + 2rikp term, which is monotonically increasing with T, 
becomes comparable in magnitude. As a result, the shear viscosity coefficient passes through 
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Figure 7 




FIG. 7: Upper panel: Simulation and theoretical results for the tracer diffusivities of a binary GC 
mixture as a function of fluid density along the isotherm T*=0.2 at the mole fraction xa=0.5. 
Lower panel: negative two-body contributions to the excess entropy of a GC mixture. 

a minimum and then grows with temperature. 

n 

As the density of GC fluid grows, its freezing temperature drops. [3| As a result, the initial 
drop of rj with T in the studied temperature range disappears, and rj increases monotonically 
with temperature both for p*=l and for p*=2. In the former case, the contributions of kinetic 
and potential terms to rj are comparable, while in the latter case the kinetic term is dominant 
for all temperatures except for the lowest one. 

Finally, in connection to our treatment of time-dependent friction of GC fluid via 
short-time ansatz, we remark that the short-time expansion coefficients for ri{t) have been 
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reported 35|] only up to the term of order f^. With this information, one would be able to 
construct an ansatz of the form given in Eq. f l23|) only for the time-dependent viscosity itself 
(or for its individual components defined in Eq. ( II ip ). but not for its first order MP. Given 
that ri[t) at high densities displays a pronounced biphasic behavior similar to that of Cy{t), 
such approach clearly would be inadequate, and we have not pursued it here. 

So far, our discussion was limited to a neat GC fluid. However, binary mixtures of GC 
particles also exhibit structural and dynamical anomalies, both in terms of density and mole 



fraction behavior. Both thermodynamic 16| and dynamic properties of these mixtures 



have been studied extensively and the following observations were made. For a binary 
mixture of particles of two different sizes, at intermediate and high densities, the tracer 
diffusivity of the larger species increases (and that of the smaller species decreases) either 
upon increasing the density of the mixture (at fixed mole fraction) or upon increasing the 
mole fraction of larger species (at a given density). As in the case of a neat GC fluid, a strong 
correlation was observed between these dynamical anomalies and structural anomalies, as 
quantified by two-body contributions to excess entropy of individual species. 

The aforementioned results for binary mixtures of GC particles were obtained via MD 
simulations. [11] In the present work, we attempt to calculate the same structural and dy- 
namical properties from a microscopic theory. We consider a two-component GC fluid 
containing species A and B. The pair interaction potential between particles has the form: 
(f)ij{r) = ejj exp[— (r/(Tjj)^], where i,j = A,B. In order to be able to perform a direct 
comparison between theory and simulation, we choose a set of parameters used in previous 



studies: 
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16[ asB = 0.665aAA, ctab = (0.5[cr^^ + c^|5])°■^ caa = £33, ^ab = OMAcaa- 
We also assume that the particles of two species have the same mass tua = m3. 

As in the case of a neat GC fluid, we compute gij{r), the radial distribution functions of 
the mixture, from the integral equation theory with HNC closure. The degree of translational 
structural order of the mixture is quantified via two-body part of the excess entropy given 



by: (ill 



S2 = '^XiS2i, (24) 

i 

where Xi is the mole fraction of species i, and S2i, which characterizes the degree of pair 



structural ordering surrounding particles of species i, is defined by:|ll[| 

^ = - ^ 2'KXjp drr"^ [gij{r) \ngij{r) - {gij{r) - 1)] (25) 
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In order to obtain the tracer diffusivities of the two species from MCT, one needs to 
construct the corresponding time- dependent friction functions, Ci{t), which, in turn, requires 
the knowledge of the matrix of dynamic and self-dynamic structure factors, Fij{k,t) and 
Fsij{k,t). However, computing FAB{k,t) from a continued fraction representation is highly 
problematic. [4S] because its zero-time second order time derivative is zero and the sign of 
the fourth derivative is oscillatory. Instead, one could obtain Fij{k, t) from a time- dependent 



density functional theory. 



43| However this approach emphasizes long-time hydrodynamic- 



like behavior of dynamic structure factors, while our results for a neat GC fluid indicate that 
transport anomalies can be explained on the basis of the short-time behavior of TCFs. In 
particular the results for D and Cy{t) of a neat GC fluid given by a short-time ansatz were for 
most part comparable in accuracy to the MCT results. Accordingly, instead of constructing 
MCT for a binary mixture, we have adopted a simpler approach by modeling (i{t) via 
Eq. (123!) . This approach requires the knowledge of the short-time expansion coefficients of 



velocity TCFs, 



't), up to the term of order t . The corresponding expressions have been 



reported earlier. [4^ 

We present our results for tracer diffusivities in dimensionless form given by: D* 



Di{mA/eAAO_ 



Simulation 
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A 



)^/^, while dimensionless density of the mixture is defined by p* = pcr^^. 



and theoretical results for D* as a function of fiuid density along the isotherm 
T* = fc^T/ 6^^=0.2 are given in the upper panel of Fig.[7]for the mole fraction x^=0.5. While 
the tracer diffusivity of the larger species displays the same anomalous density behavior as 
the self-diffusion coefficient of a neat GC fluid (i.e. passes through a minimum around 
p* ~0.4 and then increases with density), D*^ decreases monotonically with p* throughout 
the entire density range studied. As a result, the curves for the two tracer diffusivities cross 
at a certain intermediate density, beyond which the mobility of larger particles exceeds that 
of the smaller particles. Theoretical results agree well with MD data below the crossover 
density, while above this point ansatz overestimates the tracer diffusivity of the larger species 
and underestimates that of the smaller species. Still, all the trends in the density behavior 
of the two diffusivities are reproduced correctly by the theory. 

Given that the only input into the theory is the structural information, i.e. the pair dis- 
tribution functions gij{r) entering the expressions for the short-time expansion coefficients, 
one can expect a strong correlation between structural and dynamical anomalies of a GC 
binary mixture, similar to that observed for a neat fluid. Indeed, the lower panel of Fig. [7] 
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Figure 8 
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FIG. 8: Upper panel: Simulation and theoretical results for the tracer diffusivity of the larger 
species of a binary GC mixture as a function of fluid density along the isotherm T*=0.2 for three 
values of the mole fraction of A. Lower panel: negative two-body contributions to the excess entropy 
of the larger species of a binary GC mixture. 



displays the two-body contributions to the excess entropy given by Eqs. ( 12^ and fl25|) . and 
one sees that —S2a and —S2b show markedly different density behavior. In analogy to the 
diffusivity, —S2a as a function of p* follows the same pattern as —S2 of a neat fluid, namely, it 
initially grows with density, passes through a maximum, and then decreases, reflecting more 



pronounced overlaps of larger particles at higher densities. [HI] By contrast, —S2b increases 
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Figure 9 
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FIG. 9: Upper panel: Simulation and theoretical results for the tracer diffusivity of the smaller 
species of a binary GC mixture as a function of fluid density along the isotherm T*=0.2 for three 
values of the mole fraction of A. Lower panel: negative two-body contributions to the excess entropy 
of the smaller species of a binary GC mixture. 



g o f the curves for the two species. 



ll| indicating that HNC closure is 



monotonically with p, which again resuhs in the crossing 
Theoretical results for S2i agree well with the MD data, 
reliable for calculating the structure of GC mixtures. 

Having analyzed the structural and dynamical anomalies of an equimolar GC mixture, 
we now consider the effect of changing the mole fraction of larger particles on the tracer 



22 



diffusivities and structural order metrics of the two species. Tracer diffusivity of the larger 
species is shown in the upper panel of Fig. [8] as a function of fluid density along the isotherm 
T*=0.2 for several values of xa] the corresponding results for —S2a are shown in the lower 
panel. Tracer diffusivity and two-body excess entropy of the smaller species at the same 
conditions are shown in the upper and lower panels of Fig. [HI respectively. One sees that in 
the low density regime, increase in xa results in decreasing tracer diffusivity and increa sing 



structural order metric for both species. As pointed out in the earlier simulation study, [11 1 
this behavior can be expected, because increasing mole fraction of the larger species results 
in higher packing fraction of the fluid. By contrast, at higher densities, larger values of xa 
correspond to more overlaps between particles, which translates into lower values of —S2i and 
higher values of Di for both species. While the trends in the mole fraction behavior of —S2i 
and Di are the same for both species, the switch from expected to anomalous dependence 
on Xa occurs earlier (i.e. at a lower value of bulk density) for the larger species. Once again, 
theory successfully captures all the trends observed in simulations, even though there are 
minor quantitative discrepancies, mostly at high densities. 

As a final remark on GC mixtures, we note that for a given density and temperature one 
might expect anomalous (i.e. nonlinear) dependence of shear viscosity on the mole fraction 



Xa- 



43| Given that our short-time ansatz is inadequate for modeling time- dependent viscosity, 
we cannot verify this conjecture presently. One possible approach to this problem would 
be to obtain dynamic structure factors Fij{k,t) for the mixture from a fully self-consistent 



MCT framework. 45(| This will be the subject of future investigation. 



IV. CONCLUSION 

In summary, we have reported a theoretical study of structural and dynamical anoma- 
lies of a neat GC fluid and GC binary mixture. As has been discussed previously in the 



literature, l9- 



11| essentially all the observed anomalies can be traced back to the fact that 
GC potential is bounded and, as such, allows interparticle overlap, which becomes more 
and more prominent at higher densities. We have employed integral equation theory to 
compute the GC liquid structure and both MCT approach and short-time ansatz for MFs 
to obtain its transport coefficients. Theory was successful in reproducing anomalous density 
behavior of both diffusion coefficient and shear viscosity, where the former increases and the 
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latter decreases with density. One major difference between the two is that viscosity displays 
anomalous behavior over a limited density range only, following which it again increases with 
density, while diffusion anomaly persists over the entire density range studied. This fact was 
rationalized by noting that viscosity can be split into potential and kinetic contributions, 
with the former displaying anomalous density dependence and the latter behaving normally. 
With increasing density and/or temperature, the relative importance of the kinetic term 
increases, and viscosity reverts from anomalous to normal density behavior. 

A deep connection between structural and dynamical anomalies uncovered in earlier 



studies 



has been reconfirmed by noting that both density and mole fraction anomalous 
behavior of diffusion could be described via a short-time ansatz for a time-dependent friction, 
with the latter constructed from structural data only. One remaining open question concerns 
the mole fraction dependence of shear viscosity of a GC mixture, it will be the subject of 
future research. 
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